Heavy element nucleosynthesis in a collapsar 



Shin-ichirou Fujimoto 1 , Masa-aki Hashimoto 2 , Kei Kotake 3 , and Shoichi Yamada 3 ' 4 



O 
O 

(N 

X> 
PL. 

(N 



> 
O 

O 
O 

Oh; 
6 

in 



13 



ABSTRACT 

We investigate synthesis of heavy elements in a collapsar. We have calculated detailed compo- 
sition of magnetically driven jets ejected from a collapsar, which is based on long-term, magneto- 
hydrodynamic simulations of a rapidly rotating massive star of 40M Q during core collapse. We 
follow evolution of abundances of about 4000 nuclides from the collapse phase to the ejection 
phase through the jet generation phase with use of two large nuclear reaction networks. We find 
that the r-process successfully operates in the jets, so that U and Th are synthesized abundantly 
when the progenitor has large magnetic field of 10 12 G and rapidly rotating core. Abundance 
pattern inside the jets is similar compared to that of r-elements in the solar system. Heavy 
neutron-rich nuclei ~ 0.01Af Q can be ejected from the collapsar. The detailed abundances de- 
pend on nuclear properties of mass model, /3-decay rate, and fission, for nuclei near the neutron 
drip line. Furthermore, we find that p-nuclei are produced without seed nuclei: not only light p- 
nuclei, such as 74 Se, 78 Kr, 84 Sr, and 92 Mo, but also heavy p-nuclei, 113 In, 115 Sn, and 138 La, can be 
abundantly synthesized in the jets. The amounts of p-nuclei in the ejecta are much greater than 
those in core-collapse supernovae (SNe). In particular, 92 Mo, 113 In, 115 Sn, and 138 La deficient in 
the SNe, are significantly produced in the ejecta. 



Subject headings: Accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances 
supernovae: general — MHD — methods: numerical — gamma rays: bursts 
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1. Introduction 

A major fraction of elements heavier than iron 
are considered through a rapid neutron capture 
process (r-process) in an explosive event (e.g., 
Qian 2003, and references therein). Recent obser- 
vations of metal poor stars (MPSs) have revealed 
that r-elements have been produced in an early 
phase of metal enrichment of the Galaxy (e.g., 
Truran et al. 2002). Abundance profiles of r- 
elements in some MPSs have similar patterns of 
the r-elements in the solar system. The fact may 
suggest that r-process can operate in an single site, 
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or sites with a similar condition. It should be 
noted that r-elements are likely to produce in a 
site, in which a elements as well as iron group el- 
ements are not abundantly synthesized (Qian & 
Wasserburg 2003). Moreover, heavy r-elements 
(atomic number, Z > 50) suggest to have different 
origin of light r-elements (Z > 40) (Wasserburg et 
al. 1996). 

Although various sites of r-process have been 
proposed, all the proposed scenarios have some 
deficiencies: (1) Neutrino-driven winds, which 
are launched from a new-born neutron star dur- 
ing core collapse of a massive star of 10-25M Q , 
have been consider as a probable site of the r- 
process (e.g., Terasawa et al. 2002, and ref- 
erences therein). The scenario, however, re- 
quires an unlikely heavy and compact neutron 
star and the winds are possibly less abundant in 
heavy r-elements with mass number greater than 
130 (Sumiyoshi et al. 2000). (2) It is suggested 
that r-process operates inside ejecta from an O- 
Ne-Mg core, which is produced in a final stage 
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of an 8-10Mq star and is expected to explode 
promptly during the core collapse (Wheeler et al. 
1998). The r-process is shown to operate during 
the explosion of the star if the explosion energy 
is an order of 10 50 -10 51 ergs (Wanajo et al. 2003). 
However, the star reveals to explode very weakly 
with an explosion energy of 2 x 10 49 ergs and an 
ejected mass of 0.008M Q . Moreover, the ejected 
mass of the r-elements is possibly too high for 
r-elements in the Galaxy. (3) A small iron core 
of a presupernova star of ~ 11M Q is expected to 
explode promptly (e.g., Sumiyoshi et al. 2001). 
Ejecta from the core have enough low electron 
fraction, or neutron-rich, for r-process to oper- 
ate successfully. The abundance profile of the 
ejecta can reproduce that of the solar system. 
However, it is still debatable to explode the core. 

(4) r-elements can be produced abundantly in a 
neutron-rich ejecta from the merger of two neu- 
tron stars (e.g., Frciburghaus et al. 1999). The 
composition is similar to that of the solar sys- 
tem. This scenario is, however, inconsistent to 
the r-element abundances of MPSs in light of the 
Galactic chemical evolution because of the rarity 
of the neutron star mergers (Argast et al. 2004). 

(5) Recently r-process calculations in MHD jets 
have been performed for a 13M (Nishimura et al. 
2005). Though the abundance distribution is well 
reproduced compared to the solar one, their re- 
sults depend on the degree of neutronization that 
is very uncertain for multi-dimensional calcula- 
tions. 

On the other hands, any neutron capture pro- 
cesses cannot produce the p-nuclei that is 35 neu- 
tron deficient stable nuclei with mass number 
A > 74. The nuclei are considered to be synthe- 
sized by sequences of (7, n) photodisintegrations 
of s-nuclei (s-process seeds) processed during he- 
lium core burning in a massive star. The scenario 
of the synthesis of the p-nuclei via the sequential 
photodisintegrations, or p-process, is proposed in 
the oxygen/neon layers of highly evolved massive 
stars during their presupernova phase (Arnould 
1976) and during their supernova (SN) explo- 
sion (Woosley & Howard 1978). Extensive inves- 
tigation of the p-process in core-collapse SNe has 
shown that the overall abundance profile of the p- 
nuclei reproduces that in the solar system (Rayet 
et al. 1995). The scenario, however, has conspicu- 
ous shortcomings (Rayet et al. 1995); (1) some p- 



nuclei, such as 92 Mo, 94 Mo, 96 Ru, 98 Ru and 138 La, 
are underproduced compared with those in the so- 
lar system and (2) the ratio of p-nuclei to oxygen 
is smaller than that in the solar system. 

During collapse of a massive star greater than 
35-40M©, stellar core is considered to promptly 
collapse to a black hole (Woosley & Weaver 1995; 
Heger et al. 2003). When the star has sufficiently 
high angular momentum before the collapse, an 
accretion disk is formed around the hole and jets 
are shown to be launched from the inner region 
of the disk near the hole through magnetic pro- 
cesses (Proga et al. 2003; Fujimoto et al. 2006). 
Gamma-ray bursts (GRBs) are expected to be 
driven by the jets. This scenario of GRBs is called 
a collapsar model (Woosley 1993). Assisted by 
the accumulating observations implying the as- 
sociation between GRBs and the death of mas- 
sive stars (e.g., Galama et al. 1998; Hjorth et al. 
2003; Zeh, Klose, & Hartmann 2004), the col- 
lapsar model seems most promising. For accretion 
rates greater than 0.1M Q s _1 , the accretion disk is 
so dense and hot that nuclear burning is expected 
to proceed efficiently In fact, an innermost region 
of the disk related to GRBs becomes neutron-rich 
through electron capture on nuclei (Beloborodov 
2003; Pruet, Woosley & Hoffman 2003; Fujimoto 
et al. 2003b, 2004). Nucleosynthesis inside out- 
flows from the neutron-rich disk has been investi- 
gated with steady, one-dimensional models of the 
disk and the outflows (Pruet, Woosley & Hoffman 
2003; Pruet, Thompson & Hoffman 2004; Fuji- 
moto et al. 2004, 2005). Not only neutron-rich 
nuclei (Fujimoto et al. 2004) but p-nuclei (Fuji- 
moto et al. 2005) are shown to be produced inside 
the outflows. However, abundances of the outflows 
are shown to highly depend on electron fractions of 
the outflows. The electron fractions are expected 
to change during the generation phase of the out- 
flows near the base of the outflows. Nevertheless, 
the electron fraction of the outflow at the base is 
assumed to be that of the disk at the base (Pruet, 
Thompson & Hoffman 2004; Fujimoto et al. 2004, 
2005). Therefore, in order to evaluate change in 
the electron fractions of the outflows and to ex- 
amine nucleosynthesis in a collapsar, we need non 
steady, multi-dimensional simulations of a collap- 
sar from the collapsing phase of the star to the 
ejection phase of the outflows from the star. 

In the present paper, we propose jets from a 
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collapsar, or rapidly rotating massive star during 
collapse to the black hole, as a new r-process site. 
In order to investigate heavy element synthesis in 
a collapsar, we have calculated detailed composi- 
tion of magnetically-driven jets ejected from the 
collapsar, based on long-term, magnetohydrody- 
namic (MHD) simulations of a rapidly rotating, 
magnetized massive star of 4OM during core col- 
lapse, which has been performed in Fujimoto et 
al. (2006) recently. We follow evolution of the 
electron fraction and abundances of about 4000 
nuclides from the collapse phase to the ejection 
phase through the jet generation phase, with the 
aid of two large nuclear reaction networks. 

In §2 we briefly describe a numerical code for 
MHD calculation of the collapsing star, initial con- 
ditions of the star, and properties of jets from the 
star. In §3, we firstly present Lagrangian evolution 
of ejecta through the jets and then two large nu- 
clear reaction networks, in which spontaneous and 
/3-delayed fission is taken into account. It is also 
shown in §3 that r-process operates inside the jets 
and that p-nuclei that is deficient in core-collapse 
SNe can be produced inside the jets. We discuss 
effects on r-process of our assumptions in a later 
expansion phase of the ejecta and implications of 
collapsars for the Galactic chemical evolution of 
r-elements and p-elements in §4. Finally, we sum- 
marize our results in §5. 

2. MHD Calculations of a Collapsar 

We have carried out two-dimensional, Newto- 
nian MHD calculation of the collapse of a rotating 
massive star of 40M©, whose core is assumed to be 
collapsed to a black hole promptly. We present our 
numerical models and results of the collapsar, in 
particular, the production and properties of jets. 

2.1. Input Physics and Numerical Code 

The numerical code for the MHD calculation 
employed in this paper is based on the ZEUS-2D 
code (Stone & Norman 1992). We have extended 
the code to include a realistic EOS (Kotake et 
al. 2004) based on the relativistic mean field the- 
ory (Shen et al. 1998). For lower density regime 
(p < 10 5 gcm~ 3 ), where no data is available in 
the EOS table with the Shen EOS, we use another 
EOS (Blinnikov et al. 1996). We consider neutrino 
cooling processes. The total neutrino cooling rate 



is evaluated with a simplified neutrino transfer 
model based on the two-stream approximation (Di 
Matteo, Perna, & Narayan 2002), with which we 
can treat the optically thin and thick regimes on 
neutrino reaction approximately. We ignore re- 
sistive heating, whose properties are highly un- 
certain, not as in Proga et al. (2003). Spherical 
coordinates, (r, 8, <p) are used in our simulations 
and the computational domain is extended over 
50km < r < 10000km and < 9 < tt/2 and cov- 
ered with 200(r) x 24(6>) meshes. Fluid is freely ab- 
sorbed through the inner boundary of 50km, which 
mimics a surface of the black hole, whose mass is 
continuously increased by the mass of the infalling 
gas through the inner boundary. We assume the 
fluid is axisymmetric and the mirror symmetry 
on the equatorial plane. We mimic strong grav- 
ity around the black hole in terms of the pseudo- 
Newtonian potential (Paczyhsky & Wiita 1980). 

2.2. Initial conditions 

We set the initial profiles of the density, temper- 
ature and electron fraction to those of the spheri- 
cal model of a 40M Q massive star before the col- 
lapse (Hashimoto 1995). The radial and azimuthal 
velocities are set to be zero initially, and increase 
due to the collapse induced by the central hole 
and self-gravity of the star. The computational 
domain is extended from the iron core to an in- 
ner oxygen layer and encloses about 4M Q of the 
star. The boundaries of the silicon layers between 
the iron core and the oxygen layers are located at 
about 1800 km (1.88M ) and 3900km (2.4M ), 
respectively. We adopt an analytical form of the 
angular velocity fl(r) of the star before the col- 
lapse: 

as in the previous study of a collapsar (Mizuno, 
Yamada, Koide, & Shibata 2004a,b) and SNe (Ko- 
take et al. 2004). Here fi and R are parameters 
of our model. Initial magnetic field, Bq, is as- 
sumed to be uniform, parallel to the rotational axis 
of the star. We consider case with f2 = 10 rad s _1 , 
i?o = 1000km, and £?o = 10 12 G, or case with a 
strongly magnetized progenitor whose core rotates 
rapidly, because we pay attention to nucleosyn- 
thesis of heavy elements, in particular r- and p- 
elements. As shown in later (§3.3), r-process can 
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operate successfully and some p-elements be syn- 
thesized in jets from the star during core collapse. 

2.3. Properties of Jets 

We briefly describe results of our MHD calcula- 
tion (Fujimoto et al. 2006). Lagrangian evolution 
of physical quantities inside jets, which is impor- 
tant for nucleosynthesis, is shown in §3.1, in detail. 

We find that jets can be magnetically driven 
from the central region of the star along the rota- 
tional axis; After material reaches to the black hole 
with high angular momentum of ~ 10 17 cm 2 s _1 , 
a disk is formed inside a surface of weak shock, 
which is appeared near the hole due to the cen- 
trifugal force and propagates outward slowly. 
The magnetic fields, which are dominated by the 
toroidal component, are chiefly amplified due to 
the wrapping of the field inside the disk and prop- 
agate to the polar region along the inner bound- 
ary near the black hole through the Alfven wave. 
Eventually, the jets can be driven by the tanglcd- 
up magnetic fields at the polar region near the 
hole at t = 0.20 s. 

3. Nucleosynthesis in a collapsar 

We examine nucleosynthesis in a collapsar, 
based on the results of our MHD simulations. 
We firstly show Lagrangian evolution of ejecta 
through jets in detail and then proceed nucle- 
osynthesis inside the jets. 

3.1. Lagrangian evolution of ejecta through 
jets 

In order to calculate chemical composition of 
material inside jets, we need Lagrangian evolution 
of physical quantities, such as density, tempera- 
ture, and, velocity of the material. We adopt a 
tracer particle method (Nagataki et al. 1997) to 
calculate the Lagrangian evolution of the physical 
quantities from the Eulerian evolution obtained 
our MHD calculations of a collapsar. Particles are 
initially placed from the Fe core to an inner O- 
rich layer. The numbers of the particles in a layer 
are weighted to the mass in the layer. The total 
numbers of the particles are set to be 1000, with 
which we can follow ejecta through the jets appro- 
priately. We find that 59 particles can be ejected 
via the jets. 



We have performed MHD calculation until t = 
tf (tf — 0.36 s). After tf, we assume that par- 
ticles are adiabatic and expand spherically and 
freely. Therefore, velocity, position, density, and 
temperature of a particle are set to be v(t) = vo, 
r(t) = rj + v(t)(t - t f ), p(t) = p f (r f /r(t)) 3 , 
T(t) = Tf(rf/r(t)), respectively, where vq is con- 
stant in time and set to be Vf. Here Vf,rf,pf, and 
Tf are velocity, position, density, and temperature 
of each particle at tf, respectively. We briefly dis- 
cuss how the assumption on the expansion affects 
abundance change in the expansion phase of the 
ejecta, later in §4.2. 

Figure 1 shows trajectories of ejected particles 
with low Y e , where Y e is the electron fraction when 
the temperature of the particle is 9 x 10 9 K. We find 
that the particles with low Y e are initially located 
in the Fe core (< 1800km). The trajectories of 
the particles are complex due to convective motion 
near the equatorial plane (Fujimoto et al. 2006). 

Figure 2 shows the evolution of density and 
temperature. Initial density and temperature of 
a particle is set to be those of the presupernova in 
a layer where the particle is initially located. As 
the particle falls near the hole, the density and 
temperature increase. It should be emphasized 
that some particles become enough high densities 
(> 10 10 gcm~ 3 ) and high temperatures (> 10 10 K) 
for protons to capture electrons, as we shall see 
later. For the particle of Y e = 0.10, the density 
and temperature stay at high levels during t = 0.1- 
0.2 s because of convection motion near the black 
hole (bottom panel in Figure 1). The motion is 
important for decrease in the electron fraction, as 
we shall see later (Figure 3). 

3.2. Nuclear Reaction Network 

As shown in Figure 2 (bottom panel), the 
ejecta can attain to high temperature greater than 
9 x 10 9 K near the black hole, where material 
is in nuclear statistical equilibrium (NSE). The 
abundances of the material in NSE can be ex- 
pressed with simple analytical expressions (Clay- 
ton 1968; Qian 2003), which are specified by p, 
T, and the electron fraction as in Fujimoto et al. 
(2003b). The electron fraction changes through 
electron and positron capture on nuclei. It should 
be emphasized that changes in the electron frac- 
tion through neutrino interactions can be ignored 
because the material is transparent for neutrinos 
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Fig. 1. — Trajectories of ejected particles. We 
select the particles to have various values of Y e . 
The left and right panel cover the region of 2000km 
x 4000km and 1000km x 1000km, respectively. 

even in the most dense part of the computational 
domain (Fujimoto et al. 2006). 

In the relatively cool regime of T < 9 x 10 9 K, 
NSE breaks and the chemical composition is cal- 
culated with two large nuclear reaction networks, 
NETWORK A and NETWORK B, presented in 
Nishimura et al. (2005). The networks include 
about 4000 nuclei from neutron and proton up to 
Fermium, whose atomic number, Z, is 100 (see 
Table 1 in Nishimura et al. 2005). The networks 
contain reactions such as two and three body 
reactions, various decay channels, and electron- 
positron capture. The experimental masses and 
reaction rates are adopted if available. Otherwise, 
theoretical nuclear data, such as nuclear masses, 
rates of two body reactions, and /3-decays are cal- 
culated with mass formula based on FRDM in 
NETWORK A but with that based in ETFSI in 



Fig. 2. — Evolution of density (top panel) and 
temperature (bottom panel) for the particles 
shown in Figure 1. 

NETWORK B. We note that theoretical rates of 
two body reactions related to nuclei with Z > 93 
are same in the two network. 

Spontaneous and /3-delayed fission is taken into 
accounts in the both networks. Experimental 
half lives and branching ratios of the spontaneous 
fission are taken from Horiguchi, Tachibana, & 
Katakura (1996) and Kinsey et al. (1996). The- 
oretical formula of a half life (Kodama & Taka- 
hashi 1975, eq. (23)) with empirical fission bar- 
rier (Mamdouh, Pearson, Rayet & Tondeur 1999, 
2001) is adopted for nuclei of Z < 100 with- 
out experimental evaluation of the half live. We 
note that in Nishimura et al. (2005), the half 
lives through the spontaneous fission are set to 
be 10^ 20 s for all nuclei with the mass number, A, 
greater than 254. Nuclei of Z — 100 decay through 
sequences of (3~ decays and/or fission. However, in 
our networks, the nuclei decay only through fission 
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because our networks contain nuclei up to Z = 100 
and the fi~ channel of the nuclei of Z = 100 shut 
off artificially. If we calculate the half life of the 
nuclei through the fission with the theoritical for- 
mula of Kodama & Takahashi (1975), the nuclei, 
whose half life is long compared with the dura- 
tion of postprocessing calculation, cannot entirely 
decay. Therefore, for all nuclei of Z = 100, we ten- 
tatively set the half lives through the spontaneous 
fission to be 10~ 20 s. We note that same branching 
ratios of /3-delayed fission are taken from theoret- 
ical evaluation (Staudt & Klapdor-Klcingrothaus 
1992) in the two networks although the theoretical 
ratios depend on an adopted mass formula (Panov 
et. al 2005). Empirical formula (Kodama & Taka- 
hashi 1975, eq.(5)) is adopted about decay prod- 
ucts through the fission. We note that the distri- 
bution of fission yields is asymmetric. 

3.3. Abundances of an ejected particle 

Once we obtain the density and temperature 
evolution of a particle, we can follow abundance 
evolution of the particle during the infall and ejec- 
tion, through post-processing calculation using the 
nuclear reaction network presented in the previ- 
ous subsection. The ejected particles are initially 
located in the iron core (Figure 1). As the par- 
ticles infall near the black hole, the particles be- 
come hotter than 9 x 10 9 K, above which material 
in the particles is in NSE. We therefore switch to 
the NSE code from the nuclear reaction network 
to calculate abundance change of the particle. On 
the other hand, during the ejection, the temper- 
ature of the particles decreases and the particle 
becomes cooler than 9 x 10 9 K. We switch back to 
the reaction network from the NSE code for the 
post-processing calculation of the particle compo- 
sition. We perform the post-processing calcula- 
tion of nucleosynthesis inside the jets until ~ 10 10 
yr. Therefore, almost unstable nuclei has been de- 
cayed to their corresponding stable nuclei. 

Figure 3 shows evolution of the electron frac- 
tion for low Y e particles. For t < 0.2 s, the elec- 
tron fractions decrease due to electron capture on 
nuclei, mainly protons, near the black hole, where 
density and temperature of the particles are high. 
On the other hands, for t > 0.4 s, the electron 
fractions increase via (3~ -decays of neutron-rich 
nuclei in particles with Y e < 0.4. The fractions 
stay nearly constant values at 0.2 s < t < 0.4 s, 
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Fig. 3. — Evolution of the electron fraction for 
particles shown in Figure 1. 

at which the particles are abundant in neutrons 
and long lived nuclei with the neutron number 
around 50. /3-decay processes are therefore less 
important for decrease in the electron fraction at 
0.2 s < t < 0.4 s. We find that there exist parti- 
cles with Y e lower than 0.3. The particles can be 
so neutron-rich that r-process operates efficiently. 
Not only a region of high density and temperature 
but also staying in the region for a long time due 
to convective motion (Figure 1) are required for 
significant decrease in the electron fraction of the 
particles. The particle of low Y e ~ 0.1 has low 
entropy per baryon (~ 10-20 fee) where fee is the 
Boltzmann constant. The condition is similar to 
that of ejecta from a massive star of ~ 11M Q , in 
which the stellar material is expected to be ex- 
ploded promptly during SN explosion (Sumiyoshi 
ct al. 2001). The similar condition is also realized 
in a jet-like explosion from an inner core of a ro- 
tating star of 13M Q (Nishimura et al. 2005). The 
r-process can operate in the explosion (Nishimura 
et al. 2005) but the minimum Y e in the explosion 
(= 0.158) is higher than that in the present study 
(= 0.10). Therefore, U and Th cannot be synthe- 
sized abundantly in the explosion, not as in the 
present study, as we shall see later (Figure 12). 

We have adopted weak interaction rates as 
those in Fuller, Fowler, & Newman (1980, 1982). 
However, recent shell-model calculations of the 
rates (Langanke et al. 2003) for nuclei in the mass 
range A = 45-65 have resulted in substantial revi- 
sions to the old rate set of Fuller, Fowler, & New- 
man (1980, 1982). In order to estimate the effects 



6 



of the revisions on Y e of particles, we have cal- 
culated evolution of the electron fraction of the 
Y e = 0.10 particle using the revised weak inter- 
action rates. We find that the electron fraction 
calculated with the revised rates is slightly larger 
than that with the old rates. However, the differ- 
ence in Y e with the two rate sets is only up to two 
percents. This is because electrons are dominantly 
captured by protons and the rates of electron cap- 
ture on protons are comparable between the two 
sets of the rates in high temperature and high den- 
sity regimes. 



tion of the Y e ~ 0.1 particle. 




Mass Number 

Fig. 4. — Abundances of nuclei with respect to 
the mass number A for particles with low Y e . 
Filled squares, filled circles, open circles, and, 
open squares represent the abundances of particles 
with Y e = 0.10, 0.23, 0.30, and 0.38, respectively. 

Figure 4 shows abundances of nuclei with re- 
spect to A for particles with Y e = 0.10, 0.23, 0.30, 
and 0.38 using NETWORK B. For a lower Y e par- 
ticle, the particle becomes more neutron-rich in 
the NSE phase, and the abundances of heavier nu- 
clei increase. We find that the second-peak nuclei 
(A ~ 130) and the third-peak nuclei (A ~ 195) can 
be abundantly produced in particles with Y e less 
than 0.3 and 0.1, respectively. We also find that 
the particle with Y e = 0.10 is abundant in U and 
Th. The particle has enough low Y e for rapid neu- 
tron capture to proceed efficiently. We find that 
the nuclear flow on the nuclear chart proceeds near 
the neutron drip line and heavy neutron-rich nu- 
clei with Z > 95 are abundantly synthesized in the 
particle (Figure 5). Therefore, nuclear properties 
of mass model, /3-decay, and fission, for nuclei near 
the drip line are important for abundance calcula- 




t = 8.2e-01 sec 

Log Xn Log Xp Log Xa 
-0.65 -11.44 -3.16 



Fig. 5. — Mass fractions on the nuclear chart for 
a particle with Y e = 0.10 at t = 0.82 s calculated 
with NETWORK B. 




Mass Number 



Fig. 6. — Abundances of the particle with Y e 
= 0.10 calculated with NETWORK B (filled 
squares) and a network same as NETWORK B 
but ignored all fission processes (open squares). 

We have taken into accounts spontaneous and 
/3-delayed fission in our nuclear reaction networks. 
Fission is expected to have important role in the 
composition of particles with low Y e . Figure 6 
shows the abundances of the particle with Y e = 
0.10 calculated with NETWORK B and a network 
same as in NETWORK B but ignored all fission 
processes. We find that the abundances of nuclei 
with A > 250 as well as those with 60 < A < 
180 change significantly if the fission processes are 
ignored. This is due to decay through fission and 
fission yields, whose mass numbers are abundantly 
distributed in 60 < A < 180. We note that fission 
recycling cannot take place in the particle. 

Finally, in order to see dependence of particle 
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Fig. 7. — The abundances of particles with Y e 
= 0.10 and 0.23 calculated with NETWORK A 
and B presented in §3.2. Open squares, filled 
squares, open circles, and filled circles represent 
abundances of the particle with Y e = 0.10 (NET- 
WORK A), those with Y e = 0.10 (NETWORK B), 
those with Y e = 0.23 (NETWORK A), and those 
with Y e = 0.23 (NETWORK B), respectively. 

composition on a nuclear mass model, we present 
the composition of particles with Y e = 0.10 and 
0.23 using NETWORK A and B in Figure 7. We 
find that overall abundance profiles are similar for 
the two networks for the both particles. Detailed 
profiles are, however, rather different for the two 
networks. For the Y e = 0.10 particle, it is remark- 
able that heavy elements, such as Pb, Th, and 
U, with NETWORK B are more abundant than 
NETWORK A. 

The time-scale of a reaction sequence, which 
synthesizes heavy neutron-rich nuclei greater than 
A > 200, mainly depends on half lives of the 
waiting-point nuclei of r-process, whose neutron 
number, N, is 82 or 126 and half live is long 
compared to the ejection time-scale of the parti- 
cle. The half lives based on FRDM are shorter 
than those with ETFSI by an order of magnitude 
for nuclei with lower atomic number (Figure 8), 
such as 124 Mo (N = 82) and 190 Gd (N = 126), 
which are abundantly synthesized. Hence, suc- 
cessive neutron capture can proceed to synthesize 
heavy neutron-rich nuclei with A > 250 in NET- 
WORK A. Most of nuclei with A > 250 decay to 
not U and Th but lighter nuclei via fission (Figure 
6). Small fractions of the nuclei decay to U and 
Th through successive /3-decays firstly followed by 
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Fig. 8. — Half lives of neutron- rich nuclei with ./V 
= 82 and 126. The half lives are calculated with 
the two mass models base on FRDM (solid lines) 
and ETFSI (dashed lines). 

a-decays. On the other hand, neutron-rich nu- 
clei with A = 232, 235, and 238, can be abun- 
dantly produced in the ejecta with NETWORK 
B. Most of the nuclei decay to U and Th through 
/3-decays. Small fractions of neutron-rich nuclei 
with A > 250 also decay to U and Th via a similar 
sequence of /3 and a-decays as in case with NET- 
WORK A. This is because U and Th are underpro- 
duced in the ejecta with NETWORK A compared 
to those with NETWORK B. 

3.4. Integrated abundances 

We shall move on mass- weighted abundances in 
the ejecta through the jets. 
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Fig. 9. — Mass distribution with respect to Y e in 
ejecta. 
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Figure 9 shows mass distribution with respect 
to Y e in the ejecta. The ejected mass is 0.081M Q 
through the jets. We find that most particles have 
high Y e > 0.45, whose mass is O.O61M , and that 
masses of particles in which third-peak nuclei can 
be produced (Y e ~ 0.1) arc 2.4xl0~ 3 M Q . Masses 
of particles with Y e < 0.4 arc 1.8xlO" 2 M . 



. , , , Collapsar (FRDM) 
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Mass Number 



Fig. 10. — Integrated abundances in ejecta calcu- 
lated with NETWORK A. Filled squares and open 
squares with dotted line represent abundances of 
the collapsar and 153 Eu-scaled abundances of solar 
r-elements (Kappeler et al. 1989), respectively. 

Figure 10 shows mass weighted- abundances in 
the ejecta calculated with NETWORK A with re- 
spect to the mass number. The abundances of the 
solar r-elements (Kappeler et al. 1989) are also 
presented in Figure 10. We scale the abundances 
in the ejecta with the abundance of 153 Eu because 
most Eu, up to 94.2%, is produced via r-process 
in the solar system (Arlandini et al. 1999). We 
find that the abundances in the jets well reproduce 
the solar pattern of r-elements. It is clearly seen 
that nuclei around the second (A ~ 130) and third 
(A ~ 195) peaks can be synthesized in the ejecta 
via the jets. However, nuclei with 160 < A < 180 
are underproduced compared with those in the so- 
lar system because of the lack of a Y e = 0.15-0.20 
particle (Figure 9). Nuclei with A < 70 are also 
less abundant due to low Y e of the particles. 

Figure 11 shows mass weighted- abundances in 
the ejecta calculated with NETWORK A with re- 
spect to Z and the Eu (Z = 63)-scaled abun- 
dances of the solar r-elements (Burris et al. 2000). 
An appreciable amount of r-clemcnts of Z > 40 
can be produced inside the jets. The abundances 
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Fig. 11. — Integrated abundances in ejecta calcu- 
lated with NETWORK A. Filled squares and open 
squares with dotted line represent abundances of 
the collapsar and Eu(Z = 63)-scaled abundances 
of solar r-elements (Burris et al. 2000), respec- 
tively. 

of the jets well reproduce those in the solar sys- 
tem. However, the detailed abundance profile of 
the ejecta are rather different from that of the so- 
lar r-elements. Elements with 65 < Z < 73 and 
79 < Z < 81 are underproduced. U and Th are 
also underproduced inside the jets. 



Scaled Solar R°PE&em\Burnsj -a- 



9 BBa 



Os 



a* . ■ . ? a; 



Zr Ag 



°. y \ * ffe "a 

■ri) ; ; ■ 

Eu Pb 



Th 
U 

!. 



30 40 



Atomic number 



Fig. 12.— Same as Fig. 11 but with NETWORK 
B. 

If we use the another network, or NETWORK 
B, mass weighted- abundances in the ejecta slightly 
change but the overall profile is not largely differ- 
ent from those with NETWORK A (Figure 12). 
The abundances of the jets well reproduce those in 
the solar system, other than the ranges of Z < 40 
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and 70 < Z < 75, in which the jet abundances are 
underproduced compared with those in the solar 
system. It should be noted that U and Th, which 
are underproduced with NETWORK A, are com- 
parably produced as the solar system. 

3.5. Abundances of nuclei other than re- 
process elements 

We present abundances of nuclei other than r- 
elements. Note that our computational domain of 
MHD calculation extends over an oxygen rich layer 
( < 10000km ). Ejecta from an outer oxygen rich 
layer are not taken into account. Ejected masses 
are therefore under-estimated of elements lighter 
than iron, in particular a-elements. 

3.5.1. a- and iron group elements 

As shown in Figure 9, most particles of 
0.061M Q have high Y e > 0.45. Inside parti- 
cles with Y e <~ 0.5, a-rich freezeout takes place 
as in the ejecta via a core-collapse SN. 56 Ni, 
which decays to 56 Fe, is therefore abundantly 
produced inside the ejecta. We find that 56 Fe 
amounts to 0.013M Q , which is comparable to 
masses of 28 Si (O.O13M ), 32 S (9.4xlO- 3 M Q ), 
58 Ni (7.0xl0- 3 M Q ), and 60 Ni (6.4x 10- 3 M Q ). In- 
side particles with Y e slightly lower than 0.5, or 
~ 0.49, nuclei heavier than iron, such as 62 ' 64 Ni, 
64 ' 66 ' 68 Zn, and 78 Se, can be produced to be greater 
than 5 x 1O~ 4 M . 

3.5.2. p-elements 

We shall examine abundances and synthesis of 
p-nuclei in the ejecta through the jets. It is sug- 
gested that p-nuclei are produced inside an accre- 
tion disk with high mass accretion rates, which 
can be realized during the collapse of a rotating 
massive star (Fujimoto et al. 2003a). 

We find that p-nuclei are significantly synthe- 
sized in the ejecta in spite of neutron-richness of 
the ejecta. Figure 13 shows the mass fractions of 
p-nuclei abundantly produced in the ejecta. The 
abundances are calculated with NETWORK B. 
The fractions are normalized by those in the so- 
lar system (Anders & Grevesse 1989). We also 
show mass fractions of p-nuclei from core-collapse 
SNe averaged over an initial mass function (Rayet 
et al. 1995). We find that not only light p- nuclei, 
such as 74 Se, 78 Kr, 84 Sr, and 92 Mo, but also heavy 
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Fig. 13. — Mass fractions of p- nuclei abundantly 
produced in the ejecta through the jets. The 
fractions are normalized by those in the solar 
system (Anders & Grevesse 1989). IMF aver- 
aged abundances of ^-nuclei from core-collapse 
SNe (Rayet et al. 1995) are shown with open cir- 
cles. 

p-nuclei, 113 In, 115 Sn, and 138 La, are produced in 
the ejecta abundantly. The abundances of p-nuclei 
in the ejecta are much greater than those in core- 
collapse SNe (Rayet et al. 1995). It should be 
emphasized that 92 Mo, 113 In, 115 Sn, and 138 La, 
which are deficient in core-collapse SNe (Rayet et 
al. 1995), are significantly produced in the ejecta. 
We note that if we use NETWORK A, the abun- 
dances of heavy p-nuclei, 113 In, 115 Sn, and 138 La, 
decrease by a factor of 2, 8, and, 8, respectively. 

In order to see production mechanism of the 
light p-nuclei, we present abundance evolution in 
a tracer particle with Y e = 0.48 in Figures 14. 
We note that the particle corresponds to the par- 
ticle with Y e = 0.48 in Figure 1-3. As the par- 
ticle infalls near the black hole, the temperature 
raises up to 1.21 x 10 10 K (Figure 2). Heavy nuclei 
are destroyed through photodisintegration as the 
temperature increases. Abundances are rapidly 
enhanced of proton, neutron, and 4 He. As the 
temperature decreases during the ejection of the 
particle through the jets, heavy elements are re- 
assembled through neutron capture firstly, and 
then proton and a capture. We find that the light 
p-nuclei, 74 Se, 78 Kr, 84 Sr, and 92 Mo are produced 
via sequences of proton capture; For 92 Mo, abun- 
dances are enhanced through a reaction sequence, 
86 Kr(p,7) 87 Rb(p,7) 88 Sr(p,7) 89 Y(p,7) 90 Zr(p,7) 91 Nb(p, 7 ) 92 Mo. 
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Fig. 14. — Abundance evolution of a tracer par- 
ticle with Y e = 0.48. Mass fractions are shown 
of 86 Kr, 90 Zr, 91 Nb, and 92 Mo, whose abundances 
are multiplied by 10 5 . Fractions are also presented 
of proton, neutron, and 4 He. 

Note that 86 Kr, whose neutron number is equal to 
the magic number of 50, are abundantly synthe- 
sized through neutron-capture because of neutron- 
richness of the particle (Y e = 0.48). The reaction 
sequences, which produce the light p-nuclei, have 
been suggested to be realized in SNe (Howard et 
al. 1991). 

On the other hands, the heavy p-nuclei, 113 In, 
115 Sn, and 138 La, have different production mech- 
anism as the light p-nuclei. We find that the heavy 
p-nuclei can be produced only in a particle with 
Y e = 0.10 and through fission processes. In fact, 
the heavy p-nuclei cannot be synthesized if we ig- 
nore all decay channel through fission (Figure 6). 
138 La are produced through fission directly but 
113 In and 115 Sn are /3-decay products of fission 
yields, 113 Cd and 115 In (and their (3~ -decay par- 
ents), respectively, whose half lives are 9.3 x 10 15 
yr and 4.4 x 10 yr. We have performed the post- 
processing calculation of nucleosynthesis inside the 
jets until ~ 10 10 yr. If the calculation is performed 
for a much shorter time, say ~ 10 8 yr, the abun- 
dances are found to decrease by about two orders 
of magnitude. It should be noted that we have 
adopted the empirical asymmetric distribution of 
fission yield (Kodama & Takahashi 1975), which 
can be fit experimental yields but its extrapola- 
tion into experimentally unknown region is ques- 
tionable (Kodama & Takahashi 1975). The abun- 
dances of the heavy p-nuclei are therefore highly 



uncertain. 

4. Discussion 

4.1. Neutrino effects on nucleosynthesis 

We have ignored changes in the electron frac- 
tion through (anti-)neutrino capture on nuclei. 
This is because the material is transparent for neu- 
trinos even in the most dense part of the computa- 
tional domain, as shown in the previous study (Fu- 
jimoto et al. 2006). We have considered the col- 
lapse of a massive star to the black hole. Neutrino 
emission is therefore found to be smaller than that 
in cases of core collapse to a neutron star (Fuji- 
moto et al. 2006) because of lower neutrino emis- 
sion from the region near the black hole compared 
with that near the neutron star. Hence, jets from a 
collapsar are likely to be an appropriate site of the 
r-process because of low electron fraction and high 
neutron-richness due to low neutrino emission. 

4.2. Abundance change in an expansion 
phase 

As shown in Figure 5, rapid neutron capture 
takes place during a later phase, or t > 0.5 s, which 
is greater than tf = 0.36 s. Evolution of physical 
quantities in the phase is described in terms of 
that of adiabatic free expansion. In order to es- 
timate dependence of particle composition on the 
assumption at the phase, we calculate abundances 
of a particle with Y e = 0.10 for various expan- 
sion velocities in the later phase (§3.1). The ki- 
netic energy of the jets is larger than the internal 
and magnetic energies (Table 1). Therefore, the 
velocity of the jets is unlikely to change appre- 
ciably if the jets are adiabatic. Figure 15 shows 
the abundances for cases with vq = Vj, vq = 2vf, 
and vo = 0.5vf, which are represented with filled 
squares, open triangles, and, open circles, respec- 
tively. For slower particles, neutron capture pro- 
ceeds more efficiently. Nuclei with A ~ 170 are 
less abundant but those with A <~ 190 are en- 
hanced compared with those in a faster particle. 
Nuclei heavier than U and Th are abundantly pro- 
duced and eventually decay through fission. How- 
ever, the abundances are similar for three cases. 
We conclude that the particle composition weakly 
depends on the assumption on the later phase ex- 
pansion. 
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Fig. 15. — Abundances of the particle with Y e = 
0.10 for various expansion velocities in the later 
phase. Filled squares, open triangles, and, open 
circles represent abundances of the particle for 
vq = Vf, vq = 2vf, and vq — 0.5v/, respectively. 



4.3. Implications for Galactic chemical 
evolution 



4-3.1. r-elements 
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Fig. 16. — loge of the ejecta via the jets (open 
circles with dotted line) and CS31082-001 (filled 
circles with solid line). Note that loge are not 
scaled in the panel. 

Figure 16 shows loge^ 1 for element i of the 
ejecta through the jets, which are assumed to be 
uniformly mixed into inter stellar medium (ISM) 
of mass, Mism = 4 x 10 5 M Q . Note that Mt S m = 



MogeA = \og(N A /N u ) + 12 and [A/B] = log(N A /N B ) 
log(N\/N^)Q for elements A and B. 



4 x 1O 5 M is slightly larger than mass of ISM 
polluted by an SN with the explosion energy of 
10 51 ergs (Shigeyama & Tsujimoto 1998). Ob- 
served log ej in an extremely metal poor star, 
CS31082-001 ([Fe/H] = -2.9), is also presented in 
Figure 16. We find that loge^ of ISM mixed with 
the ejecta are similar to those in CS31082-001 for 
Z > 38. However, for Z < 30, loge; of the ISM are 
much less than those in CS31082-001. Especially, 
iron is poor in the ISM mixed with the ejecta com- 
pared with that in CS31082-001 ([Fe/H] - 5 of the 
ISM). 

Moreover, we find that the observed loge^ of 
CS31082-001 is well reproduced with those of the 
material in ISM for elements with Z ranging from 
10 to 92. (Figure 17), if we assumed that the 
ejecta via the jets is uniformly mixed to the ISM 
of Mism, which is pre-enriched by a hyparnova 
(HN) of 4OM . The ejected masses of nuclei from 
the HN are taken from model 40C in Maeda & 
Nomoto (2003), in which the star explodes as- 
pherically with extremely high explosion energy 
of 3.23 x 10 52 erg and ejects massive 56 Ni of 0.24 
Mq. We note that nickel is over abundant if we 
adopt masses of SN ejecta of a star, which explodes 
spherically and can ejecta massive iron. 
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Fig. 17. — Same as Fig. 16 but for the ejecta 
via the jets mixed with ejecta via a hyparnova of 
40M Q . 



4.3.2. p-elements 

The ejected mass from the collapsar is 0.081 
M Q . The mass is smaller than ejected mass of 
p-nuclei, M PPL , of 0.15 M Q from a core-collapse 
SN of 13 M (Rayet et al. 1995). However, the 
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p-nuclei, 74 Se, 78 Kr, 84 Sr, 92 Mo, 113 In, 115 Sn, and 
138 La, are produced much abundantly in the ejecta 
from the collapsar by 1-2 orders of magnitude com- 
pared with those in the ejecta from the SN (Fig- 
ure 13). The masses of the above p- nuclei from 
the collapsar, therefore, are comparable or greater 
than those from the SN. The collapsar could be 
important for the Galactic chemical evolution of 
the p-nuclei if the events are not rare compared 
with core-collapse SNe. 

Moreover, p-process in oxygen-rich layers is also 
expected to operate in the collapsar, as in p- 
process layers (PPLs) in core-collapse SNe. This 
is because s-process seeds of p-process exist in the 
layers due to weak s-process during He core burn- 
ing of the star and the jets propagate in the lay- 
ers to heat the layers to enough high tempera- 
tures for p-process to operate efficiently. There- 
fore masses of p-nuclei can be enhanced via the 
p-process, Ejected p-nuclei from the collapsar are 
the sum of those in the ejecta from oxygen- rich lay- 
ers with the composition similar to that in PPLs in 
core-collapse SNe and those in the jets from inner 
core, which abundantly contain p-nuclei deficient 
in the SNe. 

Although ejected mass of p-nuclei is possibly 
smaller than Mppl because of aspherical explosion 
in the collapsar, a ratio of the mass to ejected mass 
of oxygen can be much larger than a ratio in the 
SNe. This is because the jets from the inner core 
have small amounts of oxygen. We note that the 
ratio in the SNe is problematically much smaller 
compared with that in the solar system (Rayet et 
al. 1995). 

4.3.3. Event rates 

Finally, we briefly discuss on event rates of 
GRBs linked to the death of a massive star and 
jet-like explosion related to GRBs. A current GRB 
rate is estimated to be ~ 5 x 1CP 5 yr _1 in the Milky 
Way (Firmani et al. 2004). A GRB rate in an av- 
erage galaxy is comparable to that in the Galaxy 
and is much smaller than the rate of core-collapse 
SNe by a factor of several hundreds (Podsiadlowski 
et al. 2004) . The GRB rate is however sensitive to 
the jet opening angle (Podsiadlowski et al. 2004). 
The rate rapidly increases at red shift z ~ 1 and is 
gradually enhanced at higher red shift (Yonetoku 
et al. 2004; Murakami et al. 2005). However, an 
estimated star formation rate is constant at higher 



red shift, though the rate increases at z ~ 1 as the 
GRB rate (Giavalisco et al. 2004). Thus, the rela- 
tive GRB rate to core-collapse SNe enlarges at the 
higher z. 

Moreover, the jets in the current investigation 
are too slow (< 0.1c) and heavy (0.022M Q ) to in- 
duce GRB. Although this "failed GRB" is possibly 
hard to observe compared with normal GRBs, two 
failed GRBs are suggested to be associated with 
Type Ic SNe, SN2001em (Granot & Ramirez-Ruiz 
2004) and SN2002ap (Totani 2003). A ratio of 
failed GRBs to Type Ic SNe is estimated to be 
2/40 = 5% for nearby 40 Type Ibc SNe (Granot 
& Ramirez- Ruiz 2004) . Note that the ratio can be 
much greater than the above estimate, or 2/15 = 
13.3% if samples of Type Ibc SNe are limited to 
those with late-time (> 100 s) observations, or 15 
Type Ibc SNe (Granot & Ramirez-Ruiz 2004). Ac- 
cordingly, we can estimate a failed GRB rate to 
<~ 1.3 x 10 _4 yr _1 with an estimate of a Type Ibc 
rate of 1 x lO^yr" 1 (Podsiadlowski et al. 2004). 
Therefore, failed GRBs are possibly more frequent 
than normal GRBs. 

In brief, the jets from the failed GRBs have 
much importance to chemical evolution of galax- 
ies compared to normal GRBs because of greater 
ejection mass and frequency. In the early universe, 
or z > 1, the jets could have a significant contri- 
bution to chemical evolution of galaxies, in partic- 
ular heavy elements, which cannot be abundantly 
synthesized in core-collapse SNe. 

5. Summary and Conclusions 

We have calculated detailed composition of 
magnetically-driven jets ejected from a collapsar, 
based on long-term, magnetohydrodynamic simu- 
lations of a rapidly rotating massive star of 40M Q 
during core collapse. The magnetic field of the star 
before the collapse is set to be uniform and paral- 
lel to the rotational axis of the star. We consider 
case with the magnetic field of 10 12 G. We follow 
evolution of abundances of about 4000 nuclei from 
the collapse phase to the ejection phase through 
the jet generation phase with the aid of a NSE 
code and two large nuclear reaction networks. We 
summarize our conclusions as follows; 

1. The r-process successfully operates in the 
jets. Abundances of the jets have similar 
profile as those of r-elements in the solar 
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system, though the detailed abundances are 
different. U and Th can be abundantly syn- 
thesized in the jets. 

2. Massive neutron-rich nuclei ~ 0.01M Q can 
be ejected from the collapsar. 

3. The abundances highly depend on nuclear 
properties of mass model, /3-decay and fis- 
sion, for nuclei near the neutron drip line. 

4. Fission is important for abundances of nuclei 
with not only A > 250 but 60 < A < 180. 

5. Light p-nuclei, such as 74 Se, 78 Kr, 84 Sr, and 
92 Mo, are abundantly produced in the ejecta 
due to successive proton capture. Heavy p- 
nuclei, 113 In, 115 Sn, and 138 La, are abun- 
dantly synthesized in the jets through fis- 
sion. 

6. The abundances of p-nuclei in the ejecta 
are much greater than those in core-collapse 
SNe, in particular, 92 Mo, 113 In, 115 Sn, and 
138 La, which are deficient in the SNe, are 
significantly produced in the ejecta. 

7. If we assumed that the ejecta via the jets is 
uniformly mixed to the ISM of 4 x 10 5 M Q , 
which is prc-cnrichcd by a hyparnova of 
40M Q , logei of the material in the ISM 
well reproduces those of the extremely metal 
poor star, CS31082-001, for elements i with 
Z = 10-92. 
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